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Open problems in hot QCD 
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We try to give a comprehensive review of the main methods used in modern multi-loop calculations in finite- 
temperature field theory. While going through explicit examples, we point out similarities and differences with 
respect to the zero-temperature case, utilizing common techniques in a transparent way whenever possible. 



1. Introduction 

Perturbative computations in finite- 
temperature QCD are presently being pushed 
to the 4-loop level, resulting in a large number 
of sum-integrals over Feynman propagators to be 
evaluated. While the reduction problem can in 
principle be tackled by integration by parts (IBP) 
methods, for which - due to their prevalence in 
zero-temperature calculations - sophisticated 
algorithms and public computer programs are 
available by now, the problem of evaluating the 
resulting set of master sum-integrals still presents 
a formidable challenge. Concerning the reduction 
step, note, however, that at finite temperature 
(T) the number of master integrals at a given 
loop order is in principle unbounded, which can 
be seen already from the infinite number of mass- 
less 1-loop tadpoles [T] (see also Eq. (1^ below) , a 
situation completely different from that at T = 0, 
where finiteness can be proven rigorously [^. 

It is fair to say that the tools needed for a sys- 
tematic evaluation of multi-loop sum- integrals are 
by far not as evolved as those at zero T, where a 
number of powerful analytic and numerical meth- 
ods have been developed and made available, such 
as Mellin transforms, harmonic sums, difference 
equations or sector decomposition, to name a few. 

In contrast, the few sum- integrals that have 
been computed beyond the 2-loop level have been 
solved on a case-by-case basis (see e.g. |3I4I5J and 
references therein), mostly by carefully studying 
the integral at hand, disentangling (sub-) diver- 



gences by suitably tailored subtractions, and us- 
ing mixed numerical and analytic methods to ob- 
tain the finite terms. 

So to make progress with perturbative finite- 
temperature field theory, it would be most wel- 
come to utilize more zero-temperature machinery 
than just the IBP relations. In this note, as a first 
small step towards this goal, we intend to display 
the issues involved in evaluating a typical nontriv- 
ial sum-integral in a somewhat modern language, 
which allows to pinpoint parallels as well as key 
differences to the zero-T case. Working on a spe- 
cific 3-loop example, we will re-derive one known 
sum-integral {Si from [4], contributing to the 4- 
loop pressure of scalar theory) and, essentially by 
changing an index N in the computation, gen- 
eralize it to a new result, which will contribute 
to the matching coefficient g^ in the dimensional 
reduction framework of hot QCD [6 . 

After introducing some basic notation and 
defining a concrete one-parameter sum-integral 
that shall serve as the main vehicle to display the 
various techniques, we will exhibit the main tools 
and ideas needed to systematically dissect the in- 
tegral into divergent (but analytically tractable) 
and finite (but more difficult) pieces in Sec. |3l 
The following two sections deal with those (more 
difficult) pieces, deriving simple one-dimensional 
integral representations which are then evaluated 
numerically. While Sec. [6] is somewhat outside 
the main flow of the paper, and serves to make 
available a number of useful formulae. Sec. [7] con- 
tains the main new result. 
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2. Notation and preliminary remarks 



of massless 3-loop basketball-type sum-integrals 



Perturbative calculations in field theories at 
non-zero temperature (T) can be organized in 
large parts in exact analogy with zero-T ones, in 
particular for situations where the system is in 
thermal equilibrium. Key differences are the ad- 
ditional scale T involved, and the manifest break- 
ing of 4-dimensional Lorentz symmetry, both ef- 
fects being induced by the presence of a heat bath 
to which the system under study is coupled. As a 
consequence, the temporal direction is compact- 
ified on a circle, leading to discretized Fourier 
modes, which have to be summed over. Work- 
ing in dimensional regularization, this amounts 
to changing the familiar integral measure as 
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while, using Euclidean notation, four-momenta 
are written as Q = {qo,q) with Q^ = qo+(^ where 
go ~ 2TiTn with n ^TL the summation index for 
the bosonic case, to which we will stick through- 
out this note. We will write d-dimensional results 
using d = 3 — 2e, and often use 



(27r)'* 



As already mentioned in the introduction, the 
sum presents a major complication when com- 
pared to similar integrals at zero T . Indeed, the 
4d integrals are contained in the sum-integrals, 
as can be easily seen by expressing the sum as a 
contour integral 

T J2 ^(27rrn) = y^-F(z)4- 
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which holds for analytic functions F{z) which 
have no poles on the real axis. The first term, 
being independent of T, therefore contains the 
leading UV behavior of the respective integral. 
We will use this fact below in Eqs. © and ^. 
Let us now concentrate on the specific example 
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of which the special cases N G {1, 2, 3} occur as 
master integrals in perturbative corrections e.g. 
to the 3-loop pressure of hot QCD [3] , to the 4- 
loop pressure of scalar theories |4| and to 3-loop 
matching coefficients [6], respectively. 

Guided by the idea that this class of integrals, 
since they originate from diagrams with two ver- 
tices, require only one integration in coordinate 
space (as opposed to three), it seems desirable to 
perform the calculation in x-space whenever pos- 
sible. However, divergences (in e as d — j- 3) ob- 
struct this simple idea. The reason is that while 
the Fourier-transformed propagator (cf. Eq. p2|) 
ff) has a simple analytic form in d = 3, it is a 
messy object for general d. Hence, care has to 
be taken to perform suitable subtractions for the 
integral, and only transform to coordinate space 
in (IR- and UV-) finite integrals whose values are 
then needed at e = only. To isolate these finite 
parts requires a series of rearrangements, which 
we will now systematically construct, using meth- 
ods pioneered in |3]. 



3. Setup: subtractions in d dimensions 

We start by separating the dominating large-P 
behavior of the massless 1-loop propagator 
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where the leading term /3 = G(l,l,d-|-1) as given 
in Sec. |6] is simply a 4d massless 1-loop bubble, 
the second term contains the massless 1-loop tad- 
pole at finite T and carries a factor of two coming 
from two ways of routing the large external mo- 
mentum P through the propagators, such that 
the UV-subtracted remainder An(P) oc l/[p2]2 
as bo|,|pl>r. 

According to Eq. ©, the sum-integral Eq. ([1]) 



decomposes as 
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In order to split off potential IR divergences 
coming from 1/[Q^]^ — l/[(7^ + Qq]^ when qo = 
(its zero-mode), we multiply each of the three 
terms by the identity {Sqo + i^—Sqg)). To complete 
the IR-subtraction for the third term, we multiply 
5qg by the identity (Spg + {1 ~ Sp„)). To complete 
the UV-subtraction, we treat the (1 — 6qg) pieces 
of the first two terms as above: in the second term 
we can once again use Eq. ([2]) for ^p pi(p_Q)i , 
while in the first term, we perform the analogous 
decomposition 
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where the primed sums denote J2n ~ En=^o- 

The sum-integral that has been written out ex- 
plicitly in Eq. ([5]) can also be trivially solved by 
adding scale-free integrals that vanish in dimen- 
sional reduction, 
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where we have again identified the UV-leading 4d 
1-loop propagator /3 — G{^^,l,d+ 1) and the 
sub-leading behavior containing the 1-loop tad- 
pole sum-integral coming from only one routing 
of the large external momentum, such that the 
UV-subtracted remainder An(Q) oc l/Q^ in this 
case. 

The three terms of Eq. ([3|) can hence immedi- 
ately be rewritten (keeping the relative ordering 
of terms for clarity) as 4-1-4-1-3=11 terms as 
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where we have identified 1- and 2-loop vacuum 
sum-integrals Is and A{si, 82,83) for which ex- 
plicit analytic results are given in Sec. IHl and de- 
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One of the four sum- integrals B^~^^, namely 
containing a 3d 1-loop sub-integral where 
the zero-component of the external momentum 
plays the role of a mass, can be simplified us- 
ing integration-by-parts (IBP) identities that are 
so profitably employed in perturbative computa- 
tions at T = 0, as witnessed by numerous con- 
tributions at this conference. This essentially au- 
tomatizes the further IR subtractions that would 
otherwise have been necessary for this term since 
it still contains the zero-mode of its l/[(5^]^ prop- 
agator, see [1]. In fact, the IBP-reduction termi- 
nates with two master integrals, (one of which is 



the massive Hoop tadpole G(l,d) given in Sec-jS]), 
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The values bi{d,x) — 1 and ai((i, a;) ~ follow 
by definition, while utilizing e.g. FIRE |7j, one 
obtains the polynomials 

62(^,3;) = (d-3)(2x-l) 
a2{d,x) = d — 2 

b3{d,x) = -{d - 3){4:{d - 5)x{x - 1) + d - A) 
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that we need below - profiting from the observa- 
tion that 6Ar>i oc (d — 3), since we will only be 
interested in constant terms at d = 3. We will 
also make use later of the generic structure 
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with, in particular, 020 = 039 = 1 and 031 = —2. 

4. Treatment of finite terms at d = 3 

After the subtraction procedure as outlined 
above, all that is left to do is to perform the four 
remaining sum-integrals in Eqs. ©-([O]), which are 
defined such that they are convergent and hence 
only need to be evaluated in d = 3, dropping 
terms of 0{e). This will enable us to perform the 
discrete sums explicitly. Furthermore, as already 
mentioned in Sec. [21 it is now profitable to trans- 
form to coordinate space, as there will be fewer 
integrations. 



Writing the denominators in terms of their 
(spatial) Fourier transforms with the help of 
Eqs. ([22I) and ([M)) . the integration over spatial 
momenta is trivial. Using the values at d = 3 
of /3 = r(e)/(4^)2, ^ = -l/(647r2), I, = T^U 
and G(l, 3) — -~l/{An) from Sec.[6l we obtain the 
one-dimensional integral representations 
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where the coefficients ajy^n have been defined in 
Eq. ([T3|. we have started to use dimensionless 
variables f = 27rTr, pQ = po/{2ttT) and used the 
abbreviation A7r(f,po) for 

^e-(l'^°l+l''''~P''l-|P''l)''- (\po\ + ^) -- (14) 

as well as A7r(f, po) for the combination 

2^e-(l'^"l+l''°-P°l-|P°l)'' (|ro|r + l)- 
'•0 

-^ {pV + 3\po\r + 3) - ^ (bob + 1) . (15) 
r 3 

Using Eqs. pTl - dTOl) it is now straightforward 
to evaluate all sums, whence Eq. (J14p becomes in 
fact independent of |po|i smd in _B^ it is actually 
simplest to first sum over pq . This leaves us with 
one-dimensional integrals over products of poly- 
logarithms, logarithms and hyperbolic functions. 



5. Numerical evaluation for N = 2 and 3 

While it would be most desirable to obtain ex- 
pressions for _BjY~ for general iV, in practice we 
are forced to evaluate them at fixed iV. While 
a number (but not all) of them could be evalu- 
ated in closed form, it is perhaps simplest to treat 
them all on the same basis, i.e. in a numerical ap- 
proximation, for which we simply use the built-in 
routines of Mathematica 8 . 

For the special case A^ = 2 we get 



B^ = -jT-^ X 0.0269726622737(1) + 0{e) , 
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J X 0.0134942763002(1) + 0{e) , 
J X 0.0042655281176(1) + 0(e) , 
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4 X 0.0004627085472(1) + 0(e) . 



At A^ = 3, the numerical values are 
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X 0.0512974438185(1) + 0(e) , 
X 0.0163807421945(1) + ©(e) , 
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We now have all the ingredients at our disposal 
to obtain the results for the two special cases B2 
and JB3. Adding up all contributions according 
to Eq. (O, the final outcome is presented in the 
concluding section. 

6. Building blocks of the computation 

In order to not clutter the main line of deriva- 
tion with well-known expressions, let us here col- 
lect a few simple results that were used in the 
previous sections. These are mainly the formulae 
for analytically known (sum-) integrals, as well as 
our definitions for the spatial Fourier transforms. 

We have used the following zero-temperature 



integrals: the 1-loop massive tadpole 
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where Sahc... = Sc + Sf, + Sc + ■•• • 

When dealing with discrete sums, Zeta func- 
tions and polylogarithms enter through 



C(.)^5]n- , Li,(a;)^5]- (17) 



ri=l 



n=\ 



with Lii(a;) = — ln(l — a;), while hyperbolic func- 
tions appear via (m G Z) 
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where p{x) can be a polynomial. Note that 
Eq. (US]) is just a special case of Eq. ([TO)) . 

We have used the following sum-integrals: the 
1-loop bosonic tadpole 

^^ 1 2rc(2g-d) 

Q 
and a specific 2-loop tadpole 
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Finally, let us define the d-dimensional (inverse, 
spatial) Fourier transforms Ts as 
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At d = 3, 7^^ reduces to a modified Bessel function 
of second kind: 
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The particular cases that we need read fi{m) — 
/2(m) = 1, /o(to) = /3(™) = (1 + 1/™) and 
/_i(m) = /4(m) == (1 + 3/m + 3/m2). 

7. Conclusions 

Collecting all terms we finally re-derive 
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which coincides with the result presented in [4]. 
We also obtain the new result 

B3 = ItT'L't [1 + hie + b,2e' + 0{e')] 
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We have abbreviated some derivatives of negative 
zeta values as Z„' = ^p^, ^n" = ^(3^ and 
used the first of the Stieltjes constants defined 
via as) = 1/(5 - 1) + Er=o(l - srin/nl 

In closing, let us express our hope that a further 
generalization of the techniques presented here 
and possibly an automatization of the subtraction 
procedure will lead towards a computer- algebraic 
treatment of multi-loop sum-integrals, in order to 
greatly streamline higher-order perturbative com- 
putations in finite-temperature field theories. 
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